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One of the key science drivers for the development of the SKA is to observe the neutral hydrogen, 
Hi , in galaxies as a means to probe galaxy evolution across a range of environments over cosmic 
time. Over the past decade, much progress has been made in theoretical simulations and obser¬ 
vations of HI in galaxies. However, recent HI surveys on both single dish radio telescopes and 
interferometers, while providing detailed information on global HI properties, the dark matter 
distribution in galaxies, as well as insight into the relationship between star formation and the 
interstellar medium, have been limited to the local universe. Ongoing and upcoming HI surveys 
on SKA pathfinder instruments will extend these measurements beyond the local universe to in¬ 
termediate redshifts with long observing programmes. We present here an overview of the HI 
science which will be possible with the increased capabilities of the SKA and which will build 
upon the expected increase in knowledge of HI in and around galaxies obtained with the SKA 
pathfinder surveys. With the SKA1 the greatest improvement over our current measurements is 
the capability to image galaxies at reasonable linear resolution and good column density sensitiv¬ 
ity to much higer redshifts (0.2 < z < 1.7). So one will not only be able to increase the number 
of detections to study the evolution of the HI mass function, but also have the sensitivity and 
resolution to study inflows and outflows to and from galaxies and the kinematics of the gas within 
and around galaxies as a function of environment and cosmic time out to previously unexplored 
depths. The increased sensitivity of SKA2 will allow us to image Milky Way-size galaxies out 
to redshifts of z = 1 and will provide the data required for a comprehensive picture of the HI 
content of galaxies back to z ~ 2 when the cosmic star formation rate density was at its peak. 
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Hi and Galaxy Evolution 


1. Introduction 

How galaxies form and evolve is one of the fundamental questions in modem astrophysics. 
Over the past decade, multiwavelength observations of galaxies spanning a wide range of cosmic 
time have indicated that galaxy evolution seems to depend primarily on two parameters: a galaxy’s 
stellar mass, and the environment in which it is located. In the A-cold dark matter (CDM) hierar¬ 
chical structure formation picture, galaxies form through the successive mergers of smaller units 
(namely dark matter haloes) to form larger structures. Prior to the epoch of galaxy formation, 
baryons existed almost entirely in gaseous form, and it is through the infall of this material onto the 
filamentary structures of the cosmic web, its accretion into the deepest potential wells, and ultimate 
collapse into dense molecular clouds, that galaxies were able to start forming and producing stars. 
The ongoing influence of dark matter haloes, and complex feedback mechanisms between baryonic 
components, continue to regulate the gas content and evolution of galaxies today. In order to fully 
understand the build up of stellar mass in galaxies, we need to understand the role that neutral gas, 
the fuel for star formation, plays, along with the physical processes involved. Therefore we need to 
study the neutral gas distributions and kinematics of galaxies, in large numbers, in different envi¬ 
ronments, as a function of cosmic time to uncover the full picture of galaxy evolution and structure 
formation in the Universe. This is one of the main scientific drivers for the development of the 
SKA. 

Along with recent progress in observations, there have been many theoretical advances in both 
semi-analytic modelling and full hydrodynamical modelling of galaxy evolution and, importantly, 
the role of neutral hydrogen gas, in both its atomic (Hi) and molecular forms (see section 2.2). 
However, observations of neutral gas have lagged behind observations at other wavelengths; due to 
the intrinsic faintness of HI emission, unreasonably long observing times are required to probe to 
redshifts beyond z ~ 0.25 with existing facilities, although with the construction of SKA pathfinder 
instruments and significant upgrades to current facilities, deeper, more sensitive, observations will 
be possible in the very near future (see section 2.3). The SKA will allow us to probe HI emission 
in galaxies to previously unexplored depths and further back in cosmic history. The SKA2 will 
enable optical quality HI imaging of Milky Way size galaxies out to redshifts of z = 1 

Analysis of thousands of galaxies in the Sloan Digital Sky Survey (SDSS) has shown that they 
follow a bimodal distribution in colour as a function of stellar mass (Baldry et al. (2004)); most 
galaxies are either located in the so-called blue cloud or on the red sequence, indicating that star 
formation is either still ongoing or was quenched billions of years ago. Furthermore, the galaxy 
‘main sequence’ (star formation rate vs. stellar mass) shows a relatively smooth evolution over 
cosmic time from z~2.5toz~0 while over the same period, the star formation rate density of the 
Universe has dropped by more than an order of magnitude since its peak at z ~ 2 — 3 (Madau et al. 
(1998); Hopkins & Beacorn (2006)). Star formation also seems to have shifted from the more mas¬ 
sive galaxies to less massive galaxies at later times, a pattern known as cosmic downsizing (Cowie 
et al. (1996)). Understanding the role of the neutral gas that fuels star formation and enables the 
build up of stellar mass as well as the role of the environment in enabling or quenching these pro¬ 
cesses, is vital to disentangle the physics involved in the evolution of galaxies. To move forward, a 
number of key questions related to the role of HI in galaxies need to be answered such as: 

- What is the distribution and kinematics of the neutral (Hi) gas within and around galaxies, 
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both as a function of environment (i.e. groups/clusters vs. the field) and over cosmic time? 

- How much HI is there on average as a function of redshift? 

- How does the Mm of galaxies scale with their stellar/halo masses and other properties, 
e. g.,star formation rate, as a function of environment and redshift? 

- How important is gas accretion vs. merging in terms of building stellar mass? 

The process of star formation is intimately related to how individual galaxies evolve. On 
kpc-size scales, the Kennicutt-Schmidt law (Bigiel et al. (2008); Leroy et al. (2012); Calzetti et al. 
(2012)) describes the relationship between the star formation rate surface density and the molecular 
gas mass surface density. However, details of the conditions needed for star formation to begin at 
the scales of (molecular) clouds are as yet unclear. High resolution observations of the interstellar 
medium (ISM) in individual galaxies arc required in order to clarify the astrophysical processes 
taking place at these scales. This topic will be more fully discussed in a separate chapter. The 
Interstellar medium in Galaxies (de Blok et al. (2014)). Further discussion of the interaction of 
galaxies with their environment, the IGM and the distribution of neutral hydrogen in the cosmic 
web will be presented in the chapter The Intergalactic Medium and the Cosmic Web (Popping et al. 
(2014)). 

Observations of HI in absoiption along lines of sight to strong radio sources can provide in¬ 
formation on the neutral gas content of galaxies at higher redshifts than it is possible to reach 
directly with HI emission measurements. These arc key measurements for studies of galaxy evo¬ 
lution at higher redshifts, and will be discussed further in a separate chapter, Cool Outflows and HI 
absorbers (Morganti et al. (2014)) 

2. H i observations and theory: the status quo 

2.1 Recent observational results 

2.1.1 The Local Universe 

Recent years have ushered in a “golden age” of HI science that has seen a variety of HI 
surveys to study galaxy evolution take place. The recent single dish, large area blind HI surveys 
HIPASS (Barnes et al. (2001)) and ALFALFA (Giovanelli et al. (2005)) have shed light on the 
global HI properties of galaxies in general and provided the best estimates of the global HI mass 
function. Poor statistics at the very high Hi mass end and the very low Hi mass end (< 10 7 
Mq) limit the current results and do not resolve questions such as whether the HI mass function 
depends on the local enviroment (Zwaan et al. (2005); Springob et al. (2005)). Yet these surveys 
have revealed that there is indeed an extreme paucity of HI objects without stellar counterparts 
when the Hi mass detection limit is pushed well below 5 x 10 6 Mq, enabling the interesting 
discovery of a few local, gas-rich dwarf galaxies such as Leo T (Ryan-Weber et al. (2008)) and Leo 
P (Giovanelli et al. (2013)). In addition the surveys demonstrate the power of combining the Hi 
survey information with data from other wavebands (optical, UV, IR) to address questions relevant 
to galaxy evolution (Catinella et al. (2010); Schiminovich et al. (2010); Wong et al. (2009)). 

More detailed information about the kinematics and structure of the HI disks in individual 
galaxies has come from resolved observations using synthesis radio telescopes. Major surveys 
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carried out in the last decade are WHISP (van der Hulst et al. (2001)), THINGS (Walter et al. 
(2008)), LVHIS (Koribalski (2008)), Atlas-3D (Serra et al. (2012)), Little THINGS (Hunter et al. 
(2012)), VLA-ANGST (Ott et al. (2012)), VIVA (Chung et al. (2009)), SHIELD (Cannon et al. 
(2011)), FIGGS (Begum et al. (2008)), HALOGAS (Heald et al. (2011)) and BlueDisks (Wang 
et al. (2013)). Major results are the systematic inventory of the distribution of dark matter in 
galaxies, the prevalence of warps, tidal features, and signs of ram pressure stripping, and detailed 
relationships between star formation and the ISM including galactic fountains and the presence of 
extra-planar HI. A synopsis of results to date can be found in two recent reviews (Sancisi et al. 
(2008); van der Hulst & de Blok (2013)). 

Observations of HI in galaxies in very dense environments (a well studied example is the 
Virgo cluster (Chung et al. (2009)) demonstrate the importance of the environment. The growing 
inventory of detailed HI observations of galaxies in different environments is mixed, but suggest 
that even intermediate density environments in galaxy groups (Hess & Wilcots (2013)), and on the 
outskirts of clusters (Jaffe et al. (2013)) impact the HI content of galaxies. Dense environments in 
small volumes, such as Hickson compact groups, also exhibit HI deficiencies (Verdes-Montenegro 
et al. (2001)) with a suggestion that a reservoir of low density Hi might exist to feed the galaxies 
(Borthakur et al. (2011)). In very isolated environments, such as voids (Kreckel et al. (2012); Beygu 
et al. (2013)) there is a suggestion that accretion is an important process. Obtaining information 
for galaxies that are specifically selected to be extremely isolated will be very important (Espada 
et al. (2011)) provided detailed imaging is performed (Scott et al. (2014)). Having a good enough 
sample of HI in and around many galaxies in different enviroments to address the question of 
nature versus nurture will will be a prime goal for the SKA. 

Three surveys are currently pushing the present technology to redshifts beyond z = 0.2; the 
HIGHz survey (Catinella et al. (2008)) with Arecibo and two interferometric surveys, BUDHIES 
(Verheijen et al. (2007); Jaffe et al. (2013)) and CHILES (Fernandez et al. (2013)) which are using 
observing times approaching 1000 hours. To date, the most successful method of detecting Hi 
emission beyond z = 0.2 with current instrumentation is stacking based on position and redshift 
information from optical surveys. In the following sections we will briefly describe the current 
status of HI observations and then focus on the important next steps to be addressed by SKA1 and 
SKA2. 

2.1.2 Observations at low redshift (z < 0.4) 

Exploiting the exquisite sensitivity of the Arecibo radio telescope, the HIGHz survey provides 
a glimpse into the Hi properties of star-forming galaxies at z ~ 0.2 located in relatively isolated 
environments. The survey measured the HI content of 39 optically-selected galaxies with redshifts 
between 0.17 < z < 0.25; these arc all actively star-forming, disk-dominated systems with stellar 
and Hi masses larger than 10 10 M 0 . This sample includes not only the highest-redshift detections 
of HI emission from individual galaxies to date, but also some of the most HI -massive systems 
known. Despite being exceptionally large, the Hi reservoirs of these galaxies arc consistent with 
what is expected from their UV and optical properties.This, and the fact that the galaxies lie on the 
baryonic Tully-Fisher relation, suggests that HIGHz systems are rare, scaled-up versions of gas- 
rich disks in the local Universe. These observations provide important insight into the properties 
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of the massive, HI -rich systems that will likely dominate the next-generation HI surveys with the 
SKA and its pathfinder telescopes. 

Taking advantage of the upgraded Westerbork receiver system, a Blind, Ultra-Deep HI Envi¬ 
ronmental Survey (BUDHIES) has been earned out to image two cosmic volumes, both covering 
the redshift range z=0.164—0.224, centered on the galaxy clusters Abell 2192 at z = 0.187 and 
Abell 963 at z = 0.206. A963 is a massive, lensing, X-ray bright, Butcher-Oemler cluster with 
an unusual high fraction of blue galaxies in its center. The main scientific goal of BUDHIES is 
to characterize the gas content of these blue galaxies in relation to the gas content of galaxies in 
the cluster outskirts and the surrounding field in which this cluster is embedded. A 2192 is a less 
massive cluster that mainly serves as a control object. The clusters within their Abell radii occupy 
only ~5% of the total surveyed volume, which also contains voids and galaxy overdensities in the 
fore- and background of both clusters. 

Data were collected with the WSRT with 78 x 12 hr on A2192 and 117 x 12 hr on A963 to obtain 
a similar HI mass limit of 2x 10 9 M 0 . A total of 39 galaxies were detected in the volume containing 
A2192 and 120 galaxies in A963. First results from a pilot survey showed that none of the blue 
galaxies in A963 were detected, while stacking the HI spectra from galaxies of similar luminosity 
and color in the surrounding field showed a clear detection (Verheijen et al. (2007)). The full dataset 
is currently being analysed with a focus on the HI content of galaxies in the various substructures 
within and around the clusters (Jaffe et al. (2012, 2013)), as well as the Hi-based Tully-Fisher 
relation at these redshifts. 

The COSMOS Hi Large Extragalactic Survey (CHILES) is an ongoing 1000 hour program 
designed to take advantage of the upgraded capabilities on the Karl G. Jansky Very Large Array 
(VLA). The upgrade enables continuous redshift coverage for Hi in emission from z = 0 — 0.5, 
doubling the look-back time of existing HI measurements. A single pointing centered in the COS¬ 
MOS field, CHILES is expected to detect and resolve Milky Way-like galaxies in Hi (Mm = 
3 x 10 10 M 0 ) at the most distant redshift accessible, and the most gas-rich galaxies that we know 
of at z = 0, totalling over 300 HI -rich objects in the roughly 4 x 10 5 Mpc 3 surveyed. The chosen 
field intersects with known large-scale structure at low, intermediate, and high redshifts, and COS¬ 
MOS provides a wealth of complementary multi-wavelength photometry and optical spectroscopy 
enabling the study of galaxy evolution as a function of galaxy properties, environment, and cosmic 
time. CHILES will be the first survey to probe the HI content, growth, and morphology of galaxies 
over the same redshift range in which star formation undergoes “cosmic downsizing”. 

First results from 50 hours observing (Fernandez et al. (2013)) produced 33 direct detections 
in observations covering a continuous redshift range z = 0 — 0.193 over a 34' x 34' field of view, 
including three galaxies for which there were no previous spectroscopic redshifts. The most distant 
galaxy was detected at z=0.176 and has an HI mass of 8 x 10 9 M 0 . Stacking the Hi spectra of 
80 galaxies in a “wall” between z = 0.12 — 0.13 produced a detection with an average HI mass of 
1.8 x 10 9 M 0 . 

Recent results from HI stacking 

Hi observations have only begun to detect 21 cm line emission from individual galaxies beyond 
redshift z ~ 0.1. Although a number of galaxies will be individually detected at these and higher 
redshifts by the upcoming SKA pathfinder surveys (see Table 1), the most stringent constraints on 
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Figure 1: Recent compilation of cosmic Hi 
gas density measurements as functions of red- 
shift and look-back time. At z=0, the hexagon 
and triangle refer to direct observations from 
blind Hi surveys (Zwaan et al. (2005); Mar¬ 
tin et al. (2010)). All points above z = 0.4 
are damped Ly-a measurements from HST and 
SDSS (Prochaska et al. (2005); Rao et al. (2006); 
Noterdaeme et al. (2009, 2012)). Bridging the 
gap between the two are the estimates from 
Hi stacking of Parkes (diamonds)(Delhaize et 
al. (2013), WSRT (yellow squares) (Rhee et al. 
(2013)) and GMRT (magenta triangle) (Lah et 
al. (2007)) observations. The curve represents 
model predictions by Lagos et al. (2012). 


the HI content of galaxies across cosmic time will come from stacking of undetected sources. 

Stacking has become a common tool to constrain the statistical properties of a population of 
objects that lack individual detections in a survey, and has been applied to a variety of different 
astrophysical data, including FI I -line spectra (Chengalur et al. (2001); Lah et al. (2007); Verheijen 
et al. (2007); Lah et al. (2009); Fabello et al. (2011); Delhaize et al. (2013); Gereb et al. (2013)). 
This technique requires independent measurements of the redshifts of the galaxies, and yields an 
estimate of the average HI content of a sample of galaxies by co-adding their line emission (see 
description in Fabello et al. (2011)). 

Over the past decade, in an effort to probe the evolution of HI gas with cosmic time, Lah et 
al. (2007), who used the Giant Metrewave Radio Telescope (GMRT) to estimate the HI content 
of star-forming galaxies at z ~ 0.24, and Rhee et al. (2013), who targeted field galaxies at redshift 
z ~ 0.1 and z ~ 0.2 with the Westerbork Synthesis Radio Telescope (WSRT), showed that it is 
possible to constrain the cosmic HI gas density, L>///, of the Universe up to look-back times 2-4 
Gyr with HI stacking. Intriguingly, these studies have begun to bridge the gap between high- 
z damped Lyman-a observations and blind 21-cm surveys at z = 0 (Fig. 1). Determining Q/// 
is especially critical at these intermediate redshifts when the cosmic star formation rate density 
plummeted. 

Stacking is also extremely promising for galaxy evolution studies in general, as a means to 
investigate the connection between gas, star formation rate, and other galaxy properties across en¬ 
vironment and time. Application of the stacking technique to a sample of ~5000 optically-selected 
galaxies with HI observations from ALFALFA has successfully recovered the main scaling rela- 
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tions connecting the gas fractions of galaxies with their structural and star formation properties (Fa- 
bello et al. (2011)) and the effects of environment on gas content (Fabello et al. (2012)). Because of 
the large statistics needed for these studies, and because of the need to reach the gas-poor regime, 
this question is optimally addressed by HI stacking, rather than by individually detected galaxies. 
Environmental effects on the gas have also been investigated at much greater distances, using HI 
stacking at z ~ 0.2 (Verheijen et al. (2007)) and z = 0.37 (Lah et al. (2009)). 

Clearly, the upcoming HI surveys with the SKA and its pathfinders will allow us to extend all 
these studies to significantly higher redshifts, and to vastly larger optically-selected samples than 
is possible today. HI stacking will be a very powerful technique to exploit the next-generation HI 
surveys, well beyond the limits imposed by individual detections. 


2.2 Recent theoretical advances 

The two most widely used techniques to study galaxy formation in a cosmological context arc 
hydrodynamical simulations and semi-analytical models (SAMs). Hydrodynamical simulations 
follow the evolution of both the dark matter and baryons (gas and stars) simultaneously, while 
SAMs evolve the baryons in a pre-computed (usually by an /V-body simulation) dark matter back¬ 
ground. These approaches arc complementary: both approximate the key physical processes that 
arc fundamental to galaxy formation, such as star formation and feedback, but hydrodynamical 
simulations directly solve the equations of gas dynamics and the gravitational interactions between 
the gas, stars and dark matter, whereas SAMs must approximate them. SAMs arc computationally 
inexpensive, can model galaxy populations in cosmological volumes over cosmic time, and can 
be used to explore the parameter space of galaxy formation rapidly. In contrast, hydrodynamical 
simulations arc computationally expensive but can capture the detailed spatial, kinematical and dy¬ 
namical structure of gas, stars, and dark matter in galaxies within individual systems or in small 
volumes, and can be used to inform the choice of physically motivated recipes used in SAMs. Both 
approaches have shown substantial progress in the treatment of gas phases in and outside galaxies 
in the last ~5 years. 

Hydrodynamical simulations have witnessed important advances in the treatment of the inter¬ 
stellar medium (ISM) in galaxies and in metal enrichment of the circumgalactic medium (CGM). 
Non-cosmological simulations of galaxy formation have been used to study in detail processes 
such as H 2 formation and destruction, the Hi to H 2 transition in non-equilibrium chemistry, non¬ 
equilibrium thermal state, and variations in the strength of the radiation field (e.g. Pelupessy et al. 
2006, 2009; (Robertson & Kravtsov (2008); Gnedin et al. (2009); Dobbs et al. (2011); Shetty & Os- 
triker (2012); Glover & Clark (2012); Christensen et al. (2012)). This work has allowed for physical 
interpretation of the observed relation between the SFR and the H 2 mass surface densities observed 
in local galaxies. Another important area of progress has been in exploration of the mechanisms 
that drive turbulence in the ISM. (Dobbs et al. (2011) simulated a large spiral galaxy and showed 
that stellar feedback is important in setting the scale height of the galaxy disk through its effect on 
the velocity dispersion. Similarly, Shetty & Ostriker (2012) find that supernovae arc the primary 
source of turbulence in the dense environments of starburst galaxies. However, Bournaud et al. 
(2010) found in their simulation of a LMC-like dwarf galaxy that gravitational instabilities - which 
drive the velocity dispersion in their simulation - rather than stellar feedback set the scale height. 
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while stellar feedback is necessary to maintain large-scale turbulence, which initiates cascades to 
small scales and suppresses the formation of very dense, small gas clumps (see also Hopkins et al. 
(2013) and Shetty & Ostriker (2012)). In cosmological hydrodynamic al simulations , the transition 
from Hi to H 2 has an important effect on the predicted number of large column density absorbers, 
Nu > 7 x 10 21 cm 2 (Altay et al. (2011)). Extensive testing of the relation between galaxies and 
Hi absorbers at high redshifts has been presented by Rahmati & Schaye (2014), who reach con¬ 
clusions similar to SAMs (Lagos et al. (2011, 2014); Popping et al. (2014); Fu et al. (2012)), that 
a large fraction of the damped Ly-a absorbers are associated with galaxies with stellar masses 
M* < 10 9 M & (see also Dave et al. (2013); Bird et al. (2014); Rafieferantsoa et al. (2014)). 

In semi-analytical models of galaxy formation, significant effort has focused on improving the 
modelling of the ISM in galaxies. Whereas previous generations of models followed the evolution 
of only a single cold gas phase (Obreschkow et al. (2009b); Power et al. (2010)), current models 
partition this cold gas phase into HI and H 2 phases and follow the evolution of these phases in a 
self-consistent manner (e.g. Cook et al. (2010); Fu et al. (2010); Lagos et al. (2011); Popping et 
al. (2014)). These models assume gas in hydrostatic and chemical equilibrium and use either the 
empirical relation of (Blitz & Rosolowsky (2006)), which relates the FF/H I gas surface density 
ratio to the local hydrostatic pressure, or theoretical relations that calculate H 2 abundance as a 
result of the equilibrium between rates of destruction due to UV radiation and formation on dust 
grains (e.g. (Krumholz et al. (2009); Gnedin & Kravtsov (2011)). Models based on the empirical 
relation of Blitz & Rosolowsky appeal - to provide a better fit to the 2=0 HI mass function (Zwaan et 
al. (2005); Martin et al. (2010)), with good agreement extending down to HI masses of 10 6 M o h 2 
(Lagos et al. (2011), and also reproduce the observed clustering of Hi selected galaxies (Kim et 
al. (2012)). 

One of the key results of SAMs is shown in Fig. 1, which shows the very mild evolution of 
G/// from z = 2 to z = 0- This prediction relates to a key question in extragalactic astrophysics: 
what accounts for the dramatically different evolution of the star formation rate (SFR) and HI 
densities of the Universe. Lagos et al. (2011) and Lagos et al. (2014) argue that the steep decline 
observed in the SFR density with decreasing redshift is closely connected to the steep decline of 
the molecular gas density - because star formation is linked explicitly to the H 2 density in galaxies 
- whereas the HI density is predicted to evolve very weakly with redshift. Lagos et al. explain this 
trend as arising from a combination of decreasing gas fractions and increasing galaxy sizes with 
decreasing redshift, both of which act to reduce gas surface density and the hydrostatic pressure of 
the disk, which controls the rate at which stars form. Therefore, the SFR density evolution reflects 
the evolution of the surface density of gas in the galaxies that dominate the SFR in the universe 
as a function of cosmic time (see also Obreschkow & Rawlings (2009)). This interpretation also 
implies an increase in the molecular to dynamical mass ratios and FF/H 1 mass ratio with increasing 
redshift, which is consistent with observations (Geach et al. (2011); Tacconi et al. (2013)). 

Recent progress in the simulation and modelling of galaxy formation has been rapid, and has 
emphasized the critical role HI plays in driving galaxy formation and evolution. This work has 
also highlighted how pivotal H I surveys on the SKA and its pathfinders will be in revealing the 
full extent of this role, how they will make a crucial contribution to multi-wavelength surveys 
campaigns, and how they will provide - arguably - the most powerful tests of the predictions of 
our theories of galaxy formation and evolution. For example, simulations make clear predictions 
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BLIND Hi SURVEYS: WIDE-AREA vs. DEEP 

Main aims: To investigate: 

- Evolution of HI content (HI mass function) of galaxies vs. z and environment 

- Galaxy formation & missing satellite problem in Local Group 

- Correlation of HI properties of galaxies and with stellar/halo masses, star-formation rates, etc. 

- Physics governing cool gas distribution & evolution at low z 

- Baryonic Tully-Fisher relation vs. z 

- Cosmological parameters relating to gas-rich galaxies 

- Cosmic web 

CHILES (JVLA) [ ongoing\ 

Specifications: Single pointing in COSMOS field, for z < 0.45 

Medium-deep blind imaging 
(APERTIF) [proposed] 

Specifications: z < 0.25 over 500 deg 2 

Northern sky Hi shallow sur¬ 
vey (APERTIF) [proposed] 

Specifications: 8 > +27°, 0 < z < 0.26 

WALLABY (ASKAP) 

[approved] 

Specifications: —90° < 8 < +30°, z < 0.26 

DINGO (ASKAP) [approved] 

Specifications: 2 Phases: 

DEEP: z < 0.26 over 150 deg 2 

UDEEP:0.1 < z < 0.43 over 60 deg 2 

LADUMA (MeerKAT) 

[approved] 

Specifications: Single pointing encompassing ECDF-S, z <1.4 
over ~4 deg 2 , 

HI SURVEYS OF INDIVIDUAL GALAXIES / CLUSTERS 

Main aims: To investigate: 

- The range of conditions in ISM in local galaxies 

- The gas-cycle in nearby galaxies by probing high vs. low column density regions 

- The connection to local cosmic web 

- Gas-stripping processes in clusters 

- HI content vs. galaxy morphological type 

MHONGOOSE (MeerKAT) 

[approved] 

Specifications: 30 nearby galaxies with 10 5 < Mm < 1O 1() M 0 , 
Nhi(5g) ~ 10 19 cm -2 

Fornax Hi survey (MeerKAT) 
[approved] 

Specifications: Fornax cluster galaxies over 11 deg 2 down to 
M m ( 3a) ~ 5 x 10 5 Mq, Nhi(5<j) ~ 10 19 cm -2 


Table 1: Overview of upcoming Hi emission line surveys on SKA pathfinder facilities, APERTIF 
[Medium-deep blind imaging survey. Northern sky Hi shallow survey (Verheijen et al. (2009)) ], ASKAP 
[WALLABY and DINGO] (Duffy et al. (2012)) and MeerKAT [LADUMA (Holwerda et al. (2012)), 
MHONGOOSE (mhongoose.astron.nl) and the Fornax HI survey (Serra (2011))]. 


for the expected relation between galaxy properties (such as stellar mass and SFRs), environment 
(measured by halo mass and also clustering) and Hi content (Lagos et al. 2011, 2014, Dave et 
al. (2013); Kim et al. (2013); Rahmati & Schaye (2014)), as well as the clustering of Hi selected 
galaxies at different HI masses and cosmic times. These measurements will be partially accessible 
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with the ASKAP and MeerKAT pathfinders, but it will be with SKA1 that the most complete and 
powerful tests will be made. 

2.3 Upcoming SKA pathfinder HI surveys 

Various Hi surveys are planned for the SKA pathfinder instruments MeerKAT, ASKAP and 
APERTIF over the next 5-7 years, ranging from deep observations of particular galaxies or clusters 
to very wide-area, relatively shallow surveys to narrow, deep surveys. A summary of the science 
aims and specifications for those surveys focusing on galaxy evolution is presented in Table 1. 

This table is meant to be a summary of HI line work to be carried out in the next 5-10 years 
which will form a solid basis for further work (more detailed and sensitive studies of selected 
objects and areas and pushing knowledge of HI in and around galaxies to higher redshifts with the 
increased collecting area of SKA1 and eventually the SKA2). The SKA2 will not only be capable 
of providing optical quality (~ 1” resolution) HI images of Milky Way size galaxies at z = 1, but 
will also have the survey speed to cany out the billion galaxy survey that will chart out the large 
scale structure of the Universe from z = 0 to z = 3 and allow precision cosmology studies (Myers 
et al. (2009)). 

3. Probing galaxy evolution with the SKA 

Once SKA1 comes online, the various surveys mentioned above will have pushed our knowl¬ 
edge of H I in galaxies and its role in galaxy evolution to a situation where we have a good inventory 
of detailed knowledge of the nearby universe. Shallow surveys will have covered most of the sky 
out to z ~ 0.25 (Duffy et al. (2012)) and provided HI masses of half a million galaxies, of which 
several tens of thousands of galaxies will be reasonably well resolved. Deeper surveys will have 
pushed to z ~ 0.6, but lack the detailed information needed to characterize the physical processes 
that govern galaxy evolution from the HI point of view: gas accretion, gas removal by winds, ram 
pressure or tidal stripping. Such details are also an important ingredient for comparisons with and 
input to the advanced modelling described in section 2.2. The obvious next steps are to obtain more 
sensitive and more detailed observations in the nearby universe (0 < z < 0.25) for characterizing 
these processes and to push single-galaxy HI mass detections to larger distances and lower HI 
masses to explore the dependence on cosmic time and environment of the HI mass function and 
other scaling relations. In the following subsections, we briefly describe these science goals in 
more detail. 

3.1 The Hi mass function 

The HI mass function (Rao & Briggs (1993)) is one of the most important statistical measures 
for gas-rich galaxy populations. It describes the volume density of galaxies as a function of HI 
mass and is therefore the neutral hydrogen equivalent of the optical luminosity function. One 
important use of the HI mass function in the context of galaxy evolution is its second moment, 
which gives the cosmic Hi mass density: the number of neutral hydrogen atoms per unit volume. 
This is an important input into the understanding of the evolution of the cosmic star formation rate 
density. When combined with damped Ly a measurements at higher redshifts (e.g., Noterdaeme et 
al. 2012, and references therein), the evolution of the HI mass density can be charted and compared 
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Figure 2: Models from Lagos et al. (2012) showing evolution of the Hi mass function with redshift and 
halo mass. Halo mass bins correspond to field (left panel), group (middle panel), and cluster (right panel) 
environments. 


to models of gas and galaxy evolution. The FI I mass density can be compared with the predictions 
of semi-analytic or hydrodynamic models in order to study the processes governing the distribution 
and evolution of cool gas (e.g., Obreschkow & Rawlings (2009), Lagos et al. 2011). 

Of equal importance is the shape of the FI I mass function. For example, Kim et al. (2012) 
demonstrate that the FI I mass function is a more sensitive probe of cosmological reionization for 
galaxy formation than the optical galaxy luminosity function. For semi-analytic models, the HI 
mass function shape at z = 0 is an indispensable constraint (Obreschkow & Rawlings (2009), Lagos 
et al. 2010, Neistein et al. 2010, Kim et al. 2011). The Hi mass function puts additional limits 
on the models compared to, e.g., the stellar mass function because the correlation between halo 
mass and cold gas mass is very different from that between halo mass and stellar mass (the most 
gas-rich galaxies often have low halo mass). Furthermore, Kim et al. (2012) show how different 
star formation laws affect the shape of the HI mass function. 

The HI mass function has been well measured for the Local Universe (Zwaan et al. (2005); 
Martin et al. (2010); however, controversy still surrounds the magnitude and sign of the environ¬ 
mental density dependence of the HI mass function, possibly due to depth and cosmic variance 
issues with existing shallow surveys (Zwaan et al. (2005); Springob et al. (2005)). Upcoming pre- 
cusor surveys such as WALLABY, will detect galaxies over a ~100 times larger volume and will 
therefore be much less sensitive to cosmic variance, enabling detailed study of the variation of the 
HI mass function with environmental density. The deeper HI pathfinder surveys such as DINGO 
and LADUMA will enable measurement of the HI mass function to intermediate redshifts, z < 0.6. 

With the SKA, it will be possible to push measurements of the HI mass function to even higher 
redshifts (z ~ 1), which will in turn place even more stringent constraints on galaxy evolution 
models (see Fig. 2). It will also allow us to probe the evolution of the cosmic neutral gas density, 
Clgas, from HI emission into the redshift range where it is currently estimated using damped Ly-a 
and magnesium absorber measurements (Prochaska et al. (2005); Rao et al. (2006); Noterdaeme et 
al. (2009, 2012); Zafar et al. (2013)). 

Many questions demanding measurements of the HI masses of galaxies can also be investi- 
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gated with these observations. As functions of both environment (i.e., field vs group or cluster)(e.g., 
Hess & Wilcots (2013); Jaffe et al. (2013)) and cosmic time, we will be able to study, for lower 
HI masses than ever before, the relationships between HI mass and stellar mass, star formation 
rate, galaxy morphology, etc., as a means to uncover how galaxies evolve over time in different 
environments. 


3.2 Using HI to probe gas inflow and removal 

Galaxies form and evolve via the acquisition of gas from the intergalactic medium (IGM). 
In star-forming (blue cloud) galaxies, gas infall is expected to occur over a Hubble time to keep 
them forming stars at a relatively high rate (Fraternali & Tomassetti (2012)). Indirect evidence for 
gas infall comes from estimates of gas depletion times (Bigiel et al. (2011)), the reconstruction of 
star formation histories (Panter et al. (2007)), and chemical evolution models of the Milky Way 
(Schonrich & Binney (2009). However, at present the observational evidence for gas inflow into 
galaxies is rather scant, and it is not clear how and in what form it takes place (Sancisi et al. (2008)). 

Inflow of gas into galaxies has a long history dating back to the discovery of the high-velocity 
clouds (HVCs) in the Milky Way (Wakker & van Woerden (1997)). To date, the origin of these 
clouds remain a mystery, and their contribution to the feeding of the Galactic star formation is 
uncertain (Putman, et al. (2012)). Deep HI observations of external galaxies also show similar fea¬ 
tures in HI (Westmeier et al. (2005)). Fig. 3 shows examples of massive clouds/filaments observed 
in projection over the bright optical disk of three nearby galaxies. These filaments/complexes are 
well separated in velocity from the HI in their disks. These structures have HI masses similar to 
that of the Galactic HVC Complex C. It is unlikely that they result from a galactic fountain and 
they may be the signs of either minor mergers or the infall of cosmological filaments (Fernandez, 
Joung, & Putman (2012)). 

The typical sensitivities achieved with the current interferometers for deep HI observations 
like those of NGC2403 (Fraternali et al. (2002)) are of order 10 19 cm 2 . SKA will improve this 
sensitivity by 2 orders of magnitude, keeping the spatial resolution very high. This will have a 
tremendous impact on our understanding of gas inflow, as we will finally be able to trace the 
origin of these structures and detect gas falling into galaxies from the IGM. At present, the only 
information about low column density material around galaxies comes from absoption line studies 
towards AGNs. From the study of Lya absorbers, it is clear that there are vast reservoirs of cold 
gas surrounding both star-forming and early type galaxies (Tumlinson et al. (2013)). However their 
origin, morphology, and fate are not clear; for instance whether they are gas clouds falling into the 
galaxies or simply floating in the galactic coronae remains ambiguous. This can be explored with 
the SKA. 

Galaxies also eject gas from their disks as a consequence of stellar winds and supernova ex¬ 
plosions. The ejected gas is likely multiphase but a large fraction of it is likely cold (Melioli et al. 
(2009)) or will cool rapidly (Houck & Bregman (1990)). Observations show HI at large distance 
from the Milky Way’s disk associated with superbubble expansion (Pidopryhora et al. (2007)) and 
gas at high velocities escaping from holes in the disks of external galaxies (Boomsma et al. (2008)). 

Such HI outflows produce large layers of extra-planar gas that contain typically 10% of the HI 
mass of a galaxy, reaching up to 25% in exceptional cases (Oosterloo et al. (2007); Sancisi et al. 
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Right Ascension (J2000) R»«hl Ascension 

Figure 3: Three examples of possible accretion of Hi in nearby spirals indicated by the presence of gas 
complexes at anomalous velocities: M 101 (van der Hulst & Sancisi (1988)), NGC 2403 (Fraternali et al. 
(2002)), and M 33 (Sancisi et al. (2008)). 


(2008)). The study of the kinematics of FI I halos gives essential information about the exchange 
of material and angular momentum between galactic disks and the IGM (Marinacci et al. (2010)). 
Flowever, at the moment, the overall picture is not understood, as these low column density features 
are not seen around all galaxies and in several cases are not very extended (Heald et al. (2011)). 
Going down to a sensitivity well below xl0 18 em 2 will provide unprecented information about 
gas escaping from galaxies and how this interacts with the surrounding medium. With SKA2 we 
will eventually be able to fully trace the gas cycle from disks to the IGM and understand its crucial 
role for the evolution of star-forming galaxies. 

Gas can also be removed from galaxies via tidal interactions (Sancisi et al. (2008)) and ram 
pressure stripping in galaxy clusters (Kenney et al. (2004), see also Fig. 4). At the moment, these 
observations are limited to very few cases of mostly very massive structures (Oosterloo & van 
Gorkom (2005)) although FI I deficiencies are also found in small but dense environments, such as 
Flickson compact groups (Verdes-Montenegro et al. (2001)). SKA will change this dramatically, 
giving the possibility to study tidal interactions and minor mergers in unprecentedly large galaxy 
samples. We will be able to investigate the prevalence of these phenomena, the accretion and 
removal rates from interaction and how all these properties change as a function of environment 
and time. The study of gas removed from galaxies in dense environments will also finally solve the 
problem of star formation quenching, something that is considered of fundamental importance for 
the migration of star-forming galaxies towards the red sequence (Blanton & Moustakas (2009)). 

3.3 The TF relation 

The Tully-Fisher relation (TFr), originally describing a correlation between a galaxy’s intrin¬ 
sic luminosity and the distance-independent width of its global HI profile (Tully & Fisher (1977)), 
is the most important observed scaling relation for gas-rich, late-type galaxies. Although the TFr is 
empirical in nature, it seems astrophysically plausible that a galaxy’s luminosity, as a proxy for its 
stellar mass, should be related to its gravitational potential as traced by the rotation speed of its gas 
disk, observed through the Doppler-broadened width of its HI emission line. The presence of dark 
matter and the detailed distribution of the visible mass inside galaxies complicate this naive inter- 
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Figure 4: Three examples of galaxies undergoing gas accretion and gas removal events. From left to right: 
NGC891 (Oosterloo et al. (2007)), NGC4111 (Verheijen (2004)), and NGC4388 (Oosterloo & van Gorkom 
(2005)). The lowest contours levels are 2 x 10 20 cm 2 (red), 5 x 10 20 cm 2 (orange) and 10 x l() 20 cm 2 
(blue). It is clear that the signs of gas accretion (NGC891, NGC4111) and of gas stripping (NGC4388) 
require column density sensitivities of well below 10 2l cm 2 . 


pretation (Zwaan et al. (1995)). Nevertheless, the TFr serves as a valuable tool to obtain distances 
to galaxies in the local universe where peculiar velocities dominate over the Flubble flow (Tully 
& Courtois (2012)), while its statistical properties constrain semi-analytical and hydrodynamical 
simulations of galaxy formation and evolution. Studies of the detailed HI kinematics of regularly- 
rotating galaxies for which the extended HI gas disk is a proper tracer of the gravitational potential, 
and for which good estimates of their relative stellar masses are available, indicate that the TFr is 
fundamentally a correlation between the maximum rotational velocity of the dark matter halo and 
the embedded total baryonic mass, regardless of how this baryonic mass is distributed (McGaugh 
et al. (2000); Verheijen (2001); McGaugh (2005)). 

Major advances in photometric imaging technologies have significantly improved our mea¬ 
sures and understanding of a galaxy’s luminosity to an extent that, combined with spectroscopic 
information, it has become common practice to replace a galaxy’s luminosity in the TFr by its de¬ 
rived stellar mass, or even its baryonic mass when taking the atomic and estimated molecular gas 
masses into account (McGaugh (2005); Zaritsky et al. (2014)). Improving our understanding of the 
width of the global FI I line as a kinematic measure has received little attention so far, and even the 
most recent studies still use this quantity (Sorce et al. (2013); Zaritsky et al. (2014)). As galaxies 
are usually not spatially resolved by single-dish radio telescopes, the shape and width of the global 
FI I profile is the result of a convolution of the detailed distribution of the FI I gas in a galaxy with 
the detailed geometry and kinematics of the gas disk, including possible warps, streaming motions, 
asymmetries and the overall shapes of rotation curves, which may not always be monotonically ris¬ 
ing to a constant circular velocity (Bosma (1981a,b); Begeman (1989); Casertano & van Gorkom 
(1991); Verheijen & Sancisi (2001); Noordermeer et al. (2007); Swaters et al. (2009)). Hence, the 
width of the global HI profile is only a proxy for the rotational speed of a galaxy, while the latter 
may not even be properly defined and identified for non-flat rotation curves, and trends in rotation 
curve shapes along the TFr may introduce significant biases in scatter and slope (Verheijen (2001); 
Noordermeer & Verheijen (2007)). Imaging the geometry and kinematics of a galaxy’s extended 
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Table 2: SKA sensitivity limits (5 sigma) for HI emission 
SKA1-MID t=10h 


Resolution d (z=0.2) 

( " ) (kpc) 

d (z=0.5) 
(kpc) 

rms noise 
mJy 

dV=5 km/s 

5 sigma N(HI) 
dV=5 km/s 

5 sigma N(HI) 
dV=25 km/s 

5 sigma M(HI) 
dV=50 km/s 
z=0.02 

5 sigma M(HI) 
dV=50 km/s 
z=0.2 

5 sigma M(HI) 
dV=50 km/s 
z=0.5 

1 3.3 

6.1 

0.141 

3.92E+21 

1.96E+22 

1.20E+07 

1.56E+09 

1.32E+10 

3 9.8 

18.4 

0.118 

3.65E+20 

1.83E+21 

1.00E+07 

1.31E+09 

1.10E+10 

10 32.8 

61.2 

0.091 

2.53E+19 

1.26E+20 

7.72E+06 

1.01E+09 

8.50E+09 

30 98.3 

183.7 

0.095 

2.95E+18 

1.47E+19 

8.11E+06 

1.06E+09 

8.92E+09 


SKA1-SUR t=10h 


Resolution d (z=0.2) 

( " ) (kpc) 

d (z=0.8) 
(kpc) 

rms noise 
mJy 

dV=5 km/s 

5 sigma N(HI) 5 sigma N(HI) 
dV=5 km/s dV=25 km/s 

5 sigma M(HI) 
dV=50 km/s 
z=0.02 

5 sigma M(HI) 
dV=50 km/s 
z=0.2 

5 sigma M(HI) 
dV=50 km/s 
z=0.5 

1 3.3 

6.1 

0.551 

1.06E+22 5.28E+22 

4.67E+07 

6.11E+09 

5.14E+10 

3 9.8 

18.4 

0.361 

9.96E+20 4.98E+21 

3.06E+07 

4.00E+09 

3.37E+10 

10 32.8 

61.2 

0.323 

1.00E+20 5.01E+20 

2.74E+07 

3.58E+09 

3.02E+10 

30 98.3 

183.7 

0.588 

1.88E+19 9.38E+19 

5.00E+07 

6.53E+09 

5.50E+10 


SKA2 t=10h 


Resolution d (z=0.2) 

( " ) (kpc) 

d (z=0.8) 
(kpc) 

rms noise 
mJy 

dV=5 km/s 

5 sigma N(HI) 5 sigma N(HI) 
dV=5 km/s dV=25 km/s 

5 sigma M(HI) 
dV=50 km/s 
z=0.02 

5 sigma M(HI) 
dV=50 km/s 
z=0.2 

5 sigma M(HI) 
dV=50 km/s 
z=0.5 

1 3.3 

6.1 

0.009 

2.53E+20 1.26E+21 

7.64E+05 

9.99E+07 

8.41E+08 

3 9.8 

18.4 

0.009 

2.81E+19 1.40E+20 

7.64E+05 

9.99E+07 

8.41E+08 

10 32.8 

61.2 

0.009 

2.53E+18 1.26E+19 

7.64E+05 

9.99E+07 

8.41E+08 

30 98.3 

183.7 

0.009 

2.81E+17 1.40E+18 

7.64E+05 

9.99E+07 

8.41E+08 


HI disk with interferometers is the only observational avenue to identify the relevant kinematic 
measure for the TFr and to significantly improve our understanding of the underlying astrophysics. 

It is expected that upcoming HI imaging surveys with ASKAP, WSRT/APERTIF and MeerKAT 
will further reduce the observed scatter in the TFr by spatially resolving the HI kinematics thereby 
improving TF-based distances to galaxies in the local universe. These imaging surveys will also 
establish the statistical properties of the TFr at modest redshifts (z < 0.2) using spatially resolved 
FI I kinematics, and up to z < 0.5 using global HI line widths. It is unlikely that the SKA will 
significantly improve upon these studies. However, given the anticipated sensitivity and angular 
resolution of the SKA, it will be possible to exploit the TFr to study the evolution of galaxies over 
cosmic time by investigating the slope, scatter and zero-point as a function of redshift and cosmic 
environment, as well as global galaxy properties such as morphological type and (specific) star for¬ 
mation rate. This requires adequate angular resolution and column density sensitivity to image the 
HI kinematics in the outer gas disks of galaxies at higher redshifts, as well as survey volumes that 
encompass all cosmic environments with sufficient statistical significance to mitigate cosmic vari¬ 
ance. In order to assess the kinematical state of the largest HI disks in Mj), galaxies at a redshift of 
z = 0.5, a spatial resolution of 10 kpc would be sufficient, corresponding to an angular resolution of 
1.6 arcseconds at a frequency of 950 MHz. Based on studies of nearby galaxies, the required col¬ 
umn density sensitivity at this angular resolution would need to be ~ 5 x 10 2<) [cm -2 ]. This could 
be achievable with an ultra-deep survey using SKA 1-MID but seems out of reach for SKA-SUR. 
The volume surveyed by a single SKA1-MID pointing at z = 0.5 at 950 MHz is large enough to 
sample all cosmic environments with a sufficient number of detections to provide adequate statis¬ 
tics on the TFr. Pushing TFr studies based on resolved HI kinematics to even higher redshifts 
would require the power of SKA2, in particular given the required column density sensitivities. 
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Table 3: Strawman SKA1 surveys 



Shallow survey 

Medium deep survey 

Deep survey 

Ultra Deep Survey 

Area 

10000 deq 2 

2000 deq 2 

150 deq 2 

2 deq 2 

Integration time 

1 yr (Sur) 

1 yr (Sur) 

1 yr (Sur) 

1 yr (Mid) 

Resolution 

10" 

5" 

3" 

2" 


0.27 

0.37 

0.69 

1.7 

mean z 

0.08 

0.11 

0.19 

0.37 

nr. of detections 

2.02E+06 

1.20E+06 

4.60E+05 

4.20E+04 


4. SKA HI Strawman Surveys 

4.1 Survey scenarios and requirements 

The science goals outlined in section 3 should lead us to define the most efficient observing 
strategies for obtaining the data, given the the technical capabilities of the SKA, in the first phase 
and eventual final form. A review of the science addressed in section 3 makes clear immediately 
that it is of the greatest importance not only to have a sufficiently large number of objects for 
statistical studies probing different regions of the universe, both spatially and in cosmic time, but 
also to have the appropriate linear resolution and column density sensitivity to probe the phenomena 
of interest. For studies involving the HI mass function, it is important to push the detection limit 
to objects with small HI masses (< 10 7 M 0 ) without having to resolve the FI I distributions. For 
all other science goals it is extremely important to resolve the galaxies and be able to detect FI I 
column densities to levels well below Nm ~ 10 21 cm" 2 . This is beautifully illustrated in Fig. 4. 

Using the specifications of the SKA baseline design (Dewdney (2013)), one can derive the 
observing parameters required to reach a depth of Ahi ~ 10 20 cm 2 at reasonable resolution. Table 
2 provides basic input to these considerations. The numbers are based on the latest sensitivity 
estimates for SKA1-MID, SKA1-SUR and SKA2. Table 2 provides information for four different 
resolutions because Ajn sensitivity scales with resolution. The corresponding linear resolution is 
given for two nominal distances. In addition, Table 2 gives the HI mass limits for a 10 hour 
integration assuming a profile width of 50 km s” 1 for three nominal distances. For reference, the 
canonical value forMjjj is ~ 8 x 10 9 M 0 (Zwaan et al. (2003, 2005); Haynes et al. (2011)). 

It is clear that in order to design the right surveys for the science described in this chapter, 
trade-offs need to be made between resolution, integration time, and area covered. Table 3 provides 
an attempt to do this for SKAl. Survey areas are chosen such that the total integration time remains 
realistic, while the appropriate resolution and column density are chosen to match the science goals 
outlined in the previous sections. For example, a survey of the entire southern sky to a depth of 
Nm ~ 10 20 at a resolution of 2" would take half a century with SKAl, so a better approach is 
to cover a smaller area at lower but still useful resolution, leaving substantial improvement on a 
large Hi survey such as WALLABY to SKA2. Which surveys cover substantial gound in redshift 
space is very well illustrated in Fig. 5, showing the distribution over redshift of the four surveys 
from Table 3. The distributions represent five statistically independent realisations of a 100 deg 2 
sky model, truncated to z = 1.2, based on the S3-SAX simulation (Obreschkow et al. (2009a,b)) 
using the telescope simulations of Popping et al. (2014) which are consistent with the sensitivity 
calculations used here. In other words SKAl is the appropriate instrument for addressing follow¬ 
up science from the SKA pathfinders; SKA2 provides a major step towards pushing our observing 
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Redshift z 


Figure 5: Distribution over redshift of the four surveys given in Table 3 based on the S3-SAX 
simulation (Obreschkow et al. (2009a,b)) 


capabilities to higher redshifts, i.e., look back times. In order to obtain a good census of galaxy 
evolution and the influence of the environment over cosmic time, one has to design surveys in such 
a way that SKA1 probes the epochs 0.1 < z < 0.5 with similar linear resolution and Nm sensitivity 
as APERTIF, ASKAP and MeerKAT will probe at 0.02 < z < 0.2, while SKA2 will eventually 
provide the data required to get the full picture of HI evolution from the current epoch to z ~ 2 
when the star formation rate density in the evolving universe is at its peak. 

5. Summary and conclusions 

Our current HI view of the universe is restricted to the Local Volume and hence only informs 
us about gas-related processes in galaxies in a late state of their - cosmic evolution. Semi-analytical 
and hydrodynamical models, however, suggest significant changes in the HI properties of galaxies 
beyond the Local Universe. We are poised to start testing these models with upcoming surveys 
on the SKA pathfinder instruments over the next few years. These facilities will provide much 
improved statistics on the HI properties of galaxies and will begin to open a window on the HI 
universe beyond the Local Volume. The SKA1 (and eventually SKA2) will truly revolutionise stud¬ 
ies of HI in galaxies by enabling observations at higher angular resolutions and with sensitivities 
that are an order of magnitude greater than currently possible. Most importantly, the SKA1 will 
allow us to reach out to distances and lookback times spanning more than half the age of the uni¬ 
verse. Combined with data at other wavelengths and detailed theoretical models, HI observations 
with the SKA1 will enable us to trace the baryon cycle in galaxies from the IGM and back out 
again in different environments. The SKA2 offers the exciting prospect of directly measuring the 
HI masses of the nearest Lyman limit systems and performing precision cosmology with the billion 
galaxy survey. Ultimately, we will reveal the evolution of the distribution and kinematics of gas in 
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galaxies back to z ~ 2 when the cosmic star formation rate density was at its peak and an order of 
magnitude higher than the present, thereby building a comprehensive picture of galaxy evolution 
over cosmic time. 
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